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A model for the radial distribution function g(r-) of a square- well fluid of variable width previously 
proposed [S. B. Yuste and A. Santos, J. Chem. Phys. 101, 2355 (1994)] is revisited and simplified. 
The model provides an explicit expression for the Laplace transform of rg{r), the coefficients being 
given as explicit functions of the density, the temperature, and the interaction range. In the limits 
corresponding to hard spheres and sticky hard spheres the model reduces to the analytical solutions 
of the Percus-Yevick equation for those potentials. The results can be useful to describe in a fully 
analytical way the structural and thermodynamic behavior of colloidal suspensions modeled as hard- 
core particles with a short-range attraction. Comparison with computer simulation data shows a 
general good agreement, even for relatively wide wells. 

I. INTRODUCTION 

In the simplest model of a colloidal dispersion, the interactions among the (large) solute molecules and the excluded 
volume effects of the (small) solvent molecules, which lead to intercoUoidal solvation forces, are incorporated by 
treating the solute particles as hard spheres (HS) of diameter a. This allowS||One to take advantage of the analytical 
solution of the Percus-Yevick (PY) integral equation for the HS potentiahlilu On the other hand, it is well known 
that the solvent particles (e.g., macromolecules) can induce an effective short-range attraction between two colloidal 
particles due to an unbalanced osmotic pressure arising from depletion of the solvent particles in the region bet^geen 
the two colloidal ones. 3 □ This explains the wideapi-ead use of Baxter's model of sticky hard spheres (SHS)lI to 
represent the properties of colloidal suspensions.lail3 This interaction model represents the attractive part of the 
potential as an infinitely deep, infinitely narrow well. In addition to the hard-core diameter cr, the model incorporates 
a "stickyness" parameter (essentially the deviation of the second virial coefficient from the HS value) that can 
be understood as a measure of the temperature: the smaller the temperature the larger the stickyness. The SHS 
potential has the advantage of leading itself to an analytical solution in the PY approximationQ and in the Mean 
Spherical Approximation (MSA)B3 On the other hancL-thf; SHS potential presents two limitations. First, the system 
of monodisperse SHS is not thermodynamically stablej£il£3 however this pathology, which is not captured by Baxter's 
solution to the PY equation, can be remedied by including some degree of polydispersity in the system. The most 
important limitation of the SHS model as a representation of a realistic short-range interaction lies in the fact that 
it is unable to distinguish two situations with the same "stickyness" (i.e. same second virial coefficient) and density, 
but different temperature and/or range. 

In order to mimic the particle interaction in colloidal systepi|S|JLa a more realistic way, a simple choice is to assume 
that the particles interact via the square-well (SW) potentiaOLj 

oo, r < a 

ip{r) = — e, a < r < Xa (1) 
0, r > Act, 

where a is the diameter of the hard core, e is the well depth, and (A — l)a is the well width. The equilibrium properties 
of a SW fluid depend on the values of three dimensionless parameters: the reduced number density p* = pa^, the 
reduced temperature T* — ksT/e (fc^ being the Boltzmann constant), and the width parameter X. In the limits 
e ^ (i.e. T* oo) and/or A ^ 1, the SW fluid becomes the HS fluid. In addition, the SHS fluicH is obtained by 
taking the limits e — > oo (i.e. T* 0) and A ^ 1, while keeping the stickyness parameter — 12(A— l){e^^'^ — 1) 
constant. If we define in the original SW fluid a generalized stickyness parameter = 4(A"^ — 1) (e^^-^" — l) as 
being proportional to the deviation of the second virial coefficient from the HS value,E3 then at a given density p* the 
parameter space can be taken as the plane (t^. A). The SHS limit explores the line (tq^, A = 1) only, while the HS 
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model corresponds to the point (r^, A) = (0,1). This geometrical picture illustrates why the SW interaction model 
can be useful to uncover a much richer spectrum of values for the relevant parameters of the problem, even if the 
attraction range is relatively short. 

Despite the mathematical simplicity of the SW p ptcn tial, no analytical solution of the conventional integral equa- 
tions fcp^fluids (YBG, HNC, PY, . . . ) is knownE3l23 The mean spherical model approximation of Sharma and 
SharmaEj provides an analytical expression for the structure factor, but it is not consistent with the hard core ex- 
clusion constraint. Mnai_of the available theoretical information about the SW fluid for variable width comes from 
perturbation theory.Lj In general, perturbation theory is based on an expansion of the relevant physical quantities 
in powers of the inverse temperature. For instance, the radial distribution function g{r; p* ,T*) of the SW fluid is 
expressed as 

g{r; p*,T*) = go{r; p*) + T*-^gi{r; p*) + • • • , (2) 

where go{r\ p*) is the radial distribution function of the reference HS fluid and 51 (r; p*) represents a first order C' 
rection. Good analytical apprtwimations for g^ are known, such as Wertheim-Thiele's-apkition of the PY equation 
Verlet- Weiss parameterization,cZl or the Generalized Mean Spherical Approxiniation.E3c3 Thus, it is inJ 
gi where different versions of perturbation theory essentially differ. A few years ago. Tang and Lu (TL)E3^ 
an analytical expression (in Laplace space) for 171, based on the MSA. Comparison with Monte Carlo (MC) simula- 
tion datao showed a general good agreement for the cases considered, but the quality of the agreement worsened 
as smaller values of A and/or T* were taken. The expectation that perturbatioUj-thispry becomes less accurate as 
the well width and the temperature decrease has been already reported elsewhere. EJcj-Jtufact, perturbation theory 
tends to overestimate the critical temperature T*^ especially for small values of A — l.L^EJ All these limitations are 
significantly apparent in the SHS limit (A — > 1, T*^^ 0), in which case the expansion (H) becomes meaningless. By 
following a completely different approach, Nezbedac3 proposed to approximate rg{r) in the interval cr < r < Ac by a 
quadratic function of r, within the context of the PY theory. The coefficients of the polynomial were then determined 
analytically by imposing the continuity of the cavity function y[r) = g[r)e^^'^'^^'''^'^ and its first two derivatives at 
r = (J. NotwithstandiiTg.ihp merits of this theory, its main limitation is that it is only applicable for very narrow wells, 
typically A — 1 ^ O.Ol.Oocj In addition, Nezbeda's theory j^ils to predict a thermodynamic critical point, except at 
A = 1, in which case Baxter's solution for SHS is recovered.c3 

In order to provide a simple theory that, while including the SHS case as a special limit, could also be applicable 
to SW fluids of variable width, Yuste and Santos proposed a model by assuming an explicit functional form for the 
Laplace transform G(l) of rg{r)£^ That functional form was suggested by the exact virial expansion of the radial 
distribution functionCJ and by the property \mij.^^+ g{r) — flnite. The parameters were subsequently determined as 
functions of p* , T*, and A by imposing the condition >5'(0) — finite, where S{q) is the structure function, as well as 
the continuity of y{r) aX r — Ac. The structural properties predicted by the model showed a good agreement with 
MC simulation results, not only fcui narrow square wells, but even for relatively wide ones (A ~ 1.5) up to densities 
slightly above the critical density.a On the other hand, the model was not fully analytical because the continuity 
condition of y{r) at r = Act gives rise to a transcendent equation that must be solved nun^erieally. In fact, the exact 
solution in the case of one-dimensional SW fluids involves a similar transcendent equationO'C3 

The aim of this paper is to propose a simpler version of the model introduced by Yuste and Santos in Ref. |44| . 
While we keep the same functional structure of G{t) and enforce the conditions g{cr^) = finite and 5(0) = finite, we 
replace the transcendent equation stemming from the continuity of y{r) at r = Ac by a quadratic equation suggested 
by the SHS limit. The resulting model has therefore a degree of algebraic simplicity similar to that of the PY solution 
for SHS (but now the parameters have a A-dependence beyond the one captured by the stickyness coefficient) and 
reduces to it in the appropriate limit. The structural properties g{r) and S{q) exhibit a fairly good, agreement with 
MC simulations,L3EZl similar to that found in the original, more complicated version of the modeLo Nevertheless, it 
is in the calculation of the thermodynamic properties (which were not addressed in Ref. where the present model 
becomes especially advantageous. In particular, the isothermal compressibility is obtained as an explicit function 
of density, temperature, and well width. This allows us to get|-the A-dependence of the critical temperature T* 
and density p*. Comparison with computer simulation estimatesc3 shows that the model predictions for T*(A) are 
remarkably good, even for values of A as large as A = 1.75. On the other hand, the predicted values of /ci*(A) are 
typically 30-45% smaller than the simulation ones, a fact that can be traced back to the solution of the PY equation 
for SHS. We al^compare the predictions of the model for lie compressibility factor at A = 1.125 and A = 1.4 with 
simulation dataOcj and with the TL perturbation theory.E3 In both cases the model presents a better agreement 
than the perturbation theory, except for temperatures larger than about twice the critical temperature. 

The paper is organized as follows. The model is introduced and worked out in Sec. |l[ with some technicalities 



being relegated to the appendices. Section III deals with the comparison with simulation results and perturbation 



theory. The paper ends with a brief discussion in Sec. IV 
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II. THE MODEL 



A. Basic requirements 

The radial distribution function g(r) of a fluid is directly related to the probability of finding two particles separated 
by a distance rEj It can be measured from neutron or x-ray diffraction experiments through the static structure factor 
S{q). Both quantities are related by Fourier transforms: 



= 1 - 27rp 
where p is the number density and 



S{q) = l + p I dre-*'i--[.g(r) - 1] 
G{t) ~ Gi-t) 



(3) 

t—iq 



G(t) = / dre-''^rg{r) (4) 
Jo 

is the Laplace transform of rg{r). The isothermal compressibility of the fluid, kt = p^^ {dp/ dp)rp, is directly related 
to the long-wavelength limit of the structure function: 

XT = pkBTKT^S{Q). (5) 

Thus, all the physically relevant information about the equilibrium state of the system is contained in g{r) or, 
equivalently, in Git). 

Now we particularize to the SW interaction potential, Eq. (|l|), so g{r) = for r < a. Henceforth, we will take 
the hard-core diameter tr = 1 as the length unit and the well depth e/fcs = 1 as the temperature unit, so that the 
asterisks in p* and T* will be dropped. It is convenient to define an auxiliary function F{t) through the relation 

Git) = t 
^ ' 1 + 12T]F{t)e-t 

OC 

= Y,{~12i^r-H[F{t)re-"\ (6) 

n=l 

where rj = ^pa^ is the packing fraction. Laplace inversion of Eq. provides an useful representation of the radial 
distribution function: 

oo 

g{r) = ^(-127?)"-V„(r - n)Q{r - n), (7) 

n=l 

where is the inverse Laplace transform of and 0(r) is Heaviside's step function. Thus, the knowledge 

of F{t) is fully equivalent to that of g{r) or S{q). In particular, the value of g{r) at contact, g(l^), is given by the 
asymptotic behavior of F{t) for large t: 

g{l+)^hiO)^ \imt^F{t). (8) 

t — ^oo 

Since 5(1"'") must be finite and different from zero, we get the condition 

F{t)r^t-'^, t^oo. (9) 
On the other hand, according to Eq. (||), the behavior of G{t) for small t determines the value of 5(0): 

G(t) = t-^ + const + ^ ~ ^^^\ + o(t^). (10) 
2477 

Iggion of Eq. (jl^) into the first equality of Eq. (^) yields the first five terms in the expansion of F{t) in powers of 
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The value ol 5(0) is related to the coefficients of and by 



24tj 



24 , 
5(0) = -ry3 



d5F(t) 



dt6 



t=0. 



1 + 8r/ + 27]^ 



(11) 



(12) 



So far, all the expressions apply to any density and any hard-core potential. As is well known, in the limit of zero 
density the cavity function y{r) = g(r)e'^^^^^''^'^ is equal to 1.l2I In the special case of the SW potential this translates 
into 



lim g{r) = (1 + a;)e(r - 1) - xe{r - A), 

T/— »0 



where x = e^/^ — 1. Equation ( [l3| ) implies that 



lim F{t) = (1 + x)it-^ + t-^) - a;e-(^-i)*(At-2 + ^-3), 

?7— ►O 



(13) 



(14) 



B. Construction of the model 



Any meaningful approximation of F{t) for the SW potential must comply with Eqs. (^, (11), and (|14|). Let us 
decompose F{t) as 



F(t) = R{t) - R(t)e 



-(A-l)t 



(15) 



The model proposed by Yuste and Santo£^ consists of assuming the following rational forms for R{t) and R{t): 

Ao+Ait — Aa+Ait 



m = 



1 + 5it + 52*2 + 53*3 ' 



1 + 5it + 52*2 + 53*3 ■ 



(16) 



This forms are compatible with ([l^). In addition, condition (||) is satisfied by construction. In fact, the contact value 
is, according to Eq. (pf). 



7(1+) 



5." 



(17) 



In order to ease the proof that the HS and SHS cases are included in Eq. (^^, it is convenient to introduce the new 
parameters A = —12r]Ao, Li = —127] \Ao{\ — 1) + Ai — vli] , and L2 = —12'qAi{\ — 1). With these changes, Eqs. 
( p^ ) and (|l^) can be recast into 



1 l + A+[Li+ L2{\ - l)-i - A{\ - 1)]* -[A + L2{X - l)-i*]e-(^-i)* 



12?? 



1 + 5i* + 52*2 + 53*3 



(18) 



where we have already taken into account the property F{0) = —l/12rj, according to Eq. (|1 ![) . The model (|1S| ) 
contains six parameters to be determined. The exact expansion (11) imposes four constraints among them. Thus, we 
can express four of the parameters in terms of, for instance, A and L2. The result is 



ii = 



1 



1 + 27/ 



1 + ^77 + 277(1 + A + A2)i2 - ^??(3 + 2A + A2)A' 



(19) 



5i = 



1 + 27; 



J + 2(1 + A + A2)i2 - i(3 + 2A + A2)A' 



52 = 



2(1 + 277) 



{-1 + 7; + 2[1 - 27;A(1 + A)]L2 - [1 - r;(l + Xf]A'} 



(20) 
(21) 
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1 + 277 



(22) 



where A' = A{X — 1)^. These four parameters are hnear functions of A and L2. Taking into account Eq. (|l2|), the 
value of >S'(0) can be expressed as a quadratic function of A and L2. Its exphcit expression is given in Appendix 

Two additional constraints are still needed to determine A and £2- First note that in the zero-density limit we 
have Li 1, Si ^ 0, S2 ^ finite, 5*3 -(12r/)-^ Thus, Eq. (|l|) is consistent with Eq. (|l|) provided that 

limA = x, lim L2 = a;A(A - 1). (23) 

r|—^0 r;— »0 

In the original formulation of the model,0 the parameter A was assumed to be independent of density, so thatyl = x. 
As for L2, it was determined by imposing the (exact) continuity condition of the function y(r) at r = AJHJ'LII which 
implies 

giX-) = {l + x)giX+). (24) 

The implenientation of this condition in the model ( p^ ) leads to a transcendent equation that must be solved 
numericallyJlj In this paper we will be concerned with a simpler version of the model in which the strong condi- 
tion lh4) is replaced by a weaker one. To that end, let us introduce a parameter r as 



12 U2 



(25) 



From (|3D and (|25D it follows that 



lim r ^ [12a;A(A - 1)]-^ (26) 

The definition (|2^) is suggested by the fact that, as proved in Appendix B of Ref. Eq. ( p^ ) is equivalent to r = tshs 
in the SHS limit, namely x — > 00, A ^ 1, x{X — 1) = (12TgHs)~^ = finite. Therefore, we may expect that a simple 
prescription for r (such that r^^ — > 12a;(A — 1) in the SHS limit) can be a good substitute for the transcendent 
equation arising from (|2j), at least for wells relatively narrow. 

Since Li, S2, and £"3 are linear functions of L2, Eq. ( p5|) shows that r is the ratio of two quadratic functions of L2. 
This relation is easily inverted to get 

^ _-l + an] + a2V^ + 0^3??^ + (1 + 2??) [1 + /ji7? + (32 V'^ + PsV^ + P^ij'^] ^'"^ 

1 (27) 



1277(70 -I- 7177 -I- 7277^) 

where the expressions of the coefficients a^, and 7^ as functions of A, A, and r are given in Appendix 

So far, we are free to fix r as a function of 77, x, and A by following any criterion we wish. In particular, we can 
enforce Eq. (p^), as done in Ref. Analogously, the parameter ^(77, x, A) (which was taken A = a; in Ref. 0) can be 
freely chosen. In the simplified version of the model we consider here, we assume that both A and r are independent 
of density, so they take values needed to satisfy the requirement ([l4| ). Those values are simply 

A^x, r = [12a;A(A- 1)]"^ (28) 

Note that this t is only slightly different from the parameter rg-yv introduced below Eq. (0), both of them becoming 
identical to rgns in the limit A ^ 1. The choice (ES) closes the construction of the model. The pair distribution 
function in Laplace space is given by Eqs. (^) and (|l^, where the expressions for the coefficients are (pj|)-(p2|), (|2^), 
and (p8|). The model provides the quantity G{t) as an explicit function of the Laplace variable t, the packing fraction 
77, the well width A, and the temperature parameter x = e^/^ — 1. Since the poles of F{t) are the roots of a cubic 
equation, the inverse Laplace transforms oit[F(t)]^ are analytically derived and then the radial distribution function 
is readily obtained from the representation (|^). From Eq. ( p8| ) it follows that the relationship between the temperature 
and the parameter r is 

T = l/ln[l + T-Vl2A(A- 1)]. (29) 

The predictions of our model to first order in the packing fraction are compared with the exact results in Appendix 
As an illustrative example, Fig. |l| shows the quantity 
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A(r) lim - ^ — ^ 



(30) 



for A = 1.1 and A — 1.125 and for the temperatures T = 0.5, 0.67, and 1. The functfon A(r) is different from zero in 
the interval 1 < r < A only, where our model slightly overestimates the value of g{r). Note that the relative deviation 
of our g{r) from the exact radial distribution function is A(r)77 for small densities. Of course, the linear growth of 
this deviation with increasing density is valid in this low-density regime only, as comparison with simulation values 
at finite densities shows [cf. Figs. ^ and |^. 




1.20 



FIG. 1. Function A(r) defined by Eq. 
(dashed lines), and 1 (dotted lines). 



j) for A = 1.1 and A = 1.125 and for the temperatures T = 0.5 (solid lines), 0.67 



Appendix ^ shows that the model reduces to Wertheim-Thiele's and Baxter's analytical solutions of the PY equation 
in the HS and SHS limits, respectively. If we consider the parameter space A-r~^, then the HS potential corresponds 
to the line t^^ — (in which case the physical properties are independent of A), while the SHS potential corresponds 
to the line A = 1 (the physical properties being r-dependent). What our model does is to extend the above picture to 
the entire plane > 0, A > 1, without compromising the mathematical simplicity present in the analytical solutions 
of the PY integral equation for HS and SHS. 



III. COMPARISON WITH SIMULATION AND OTHER THEORIES 



A. Structural properties 



The structural properties obtained from the original version of the model [i.e. with L2 detepiniajed by solving the 
transcendent equation stemming from Eq. (p^j were profusely compared with simulation datau'cd'cZl in Ref. ^4 - We 
have checked that the simplified version of the model presented here [cf. (|2^), (p8|)] gives results very close to those of 
the original model, especially for narrow wells. Consequently, we will present only a brief comparison with simulations 
in this Subsection. 
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The structure of a fluid is usually determined by neutron or x-ray scattering experiments, which measure the 
structure factor S{q). This quantity is directly related in the Laplace transform G(t) by Eq. (||) and so can be 
obtained explicitly in our model. In 1984, Huang et alEU performed MC simulations of SW fluids with A = 1.1 
to reproduce the main features of the structure factor of micellar suspensions. Figure ^ shows S{q) obtained from 
simulationcZl for A = 1.1, T^^ = 0.92, and rj = 0.07, as compared with our model, the PY solutipa-iir SHS (with 
the conventional choice of as stickyness parameter), the Tang-Lu (TL) perturbation theoryj230 and the PY 
solution for hard spheres. The deviations of the simulation data from the PY-HS curve are essentially a measure 
of effects associated with the non-zero values of the square well width and of the inverse temperature. These are 
qualitatively described by the TL perturbation theory, except for small wavenumbers. This region is well represented 
by the PY-SHS curve and by our model, but the latter is better near the first maximum. The value of the critical 
temperature for A — 1.1 can be estimated to be Tc « 0.5 [cf. Table Consequently, the case considered in Fig. ^ 
corresponds to a rather hot gas. That is why the simulation data are not too far from the HS values (except in the 
region q ^ 3). In order to highlight the differences that can be expected at a smaller temperature, the case T — 0.5 is 
considered in Fig. ^. Now, the structure factor predicted by our model and by the PY-SHS solution is very different 
from that predicted by the TL perturbation theory, the latter being very close to the PY-HS solution. This is not 
surprising, since in any perturbation theory the quantities are expanded in powers of and obviously the value 
= 2 is beyond its range of applicability. On the other hand, the curves corresponding to the PY-SHS solution 
and our model are relatively close (except for a slight phase shift), as expected from the fact that the well width 
is rather small. Moreover, since the stickyness parameter has been chosen as to reproduce the correct second virial 
coefficient, the PY-SHS solution does a general good job at this very low density. 
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FIG. 3. Structure factor, S{q), corresponding to a SW fluid with A = 1.1, rj = 0.07, and T = 0.5. The hues are the resuh 

predicted by the present model ( — ), the PY equation for SHS {■■■), the TL perturbation theory ( ), and the PY equation 

for HS (— -). 

Now we consider the radial distribution function itself. As a representative exarnple of a width not extremely small, 
we take the case A = 1.125, for which MC simulations are availableo Figure |j shows g{r) for A = 1.125 at the 
thermodynamic state T = 1, p = 0.8. As seen in the figure, the original model of Ref. exhibits a remarkable 
agreement with the simulation data. We also observe that our present model captures reasonably well the behavior of 
g{r) at this high density, while the TL approximation predicts a too small contact value g{V^). The contact values are 
plotted in Fig. || as a function of 1/T for A = 1.125 and p — 0.4, 0.6, and 0.8. At 1/T = the system corresponds to 
HS and then our model and the TL theory reduce to Wertheim-Thiele's solution of the PY equation, which tends to 
underestimate 5(1"'") at high densities. As the inverse temperature increases, the TL theory predicts a linear increase 
of .g(l^) [cf. Eq. (||)] that clearly deviates from the MC data. On the other hand, the temperature dependence of 
(7(1+) is well described by both the original and the simplified versions of the model, except at 1/T = 2. 
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FIG. 4. Radial distribution function, g{r), corresponding to a SW fluid with A = 1.125, p — 0.8, and T — 1. The circles 
represent MC data,E3 the dotted line is the result obtained from the model of Ref. the solid line is the result given by the 
present model, and the dashed line is the prediction from the TL perturbation theory. 




FIG. 5. Plot of g{l+) as a function of l/T for A = 1.125 and three densities: p = 0.4 (o), p = 0.6 (□), and p = 0.8 (O). The 
symbols represent MC data taken from Ref. ^ (1/r 7^ 0) and from Ref. ^ (1/r — 0). The the dotted lines are the results 
obtained from the model of Ref. the solid lines are the results given by the present model, and the dashed lines are the 
predictions from the TL perturbation theory. 
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B. Thermodynamic properties 



The compressibility equation of state is obtained from Eqs. (^) and (Al). This routeJs preferable to the virial route 
because the latter is known to yield an unphysical critical behavior in the SHS limit.E3 From Eq. we have 



P 



pksT 



= V 



(31) 



where the density dependence of xt is given in our model by Eqs 
explicitly, it does not allow us to perform the integration in Eq. ( 
is obtained by numerical integration. 




l|) and (27). Although this dependence is known 



Ih analytically, so that the compressibility factor Z 



Before delving into the thermodynamic predictions of our model, let us compare it with the original on E in the case 
of a moderately narrow well, namely A = 1.125. Figure]^ shows the density dependence of the inverse susceptibility for 
T = 0.5, 0.67, and 0.9. As the temperature increases, our simplified model is seen to overestimate the compressibility 
of the fluid, especially for large densities. On the other hand, at T = 0.5 (which is practically the critical temperature, 
cf. Table |), both versions of the model yield undistinguishable results. This is quite encouraging, since it is in the 
domain of low temperatures (or, equivalently, high stickyness) where our model is expected to correct the deficiencies 
of perturbation theories and become useful. 




FIG. 6. Density dependence of the inverse isothermal susceptibility for A = 1.125 and three temperatures. The dashed lines 
are the results obtained from the model of Ref. W4 and the solid lines are the results from the present model 



One of the advantages of our simplified model is that it allows us to derive explicit expressions for the coordinates 
of the critical point. The spinodal line is the locus of points where the isothermal compressibility diverges. Equation 



( |Al| ) clearly shows that xt = <S'(0) ^ oo if and only if L2 ^ oo. Thus, according to Eq. (^71), the spinodal hue is 
given by a solution to the quadratic equation 70 + jiJj + 72??^ = 0. This equation has two real solutions, ?7±(t), only 
if T is a smaller than a critical value Tc given by 

1 3 + A + 2^/2A 

9-2A + A2 ■ (^2) 
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The critical temperature Tc is obtained by setting r = Tc in Eq. (p9|). If t < Tc, the smallest root, ri-{T), does not 
define the vapor branch of the spinodal line because it is also a root of the numerator of Eq. ( p7|) and so L2 remains 
finite at 77 = rj^{T). Therefore, only the liquid branch of the spinodal line, ?y+(T), exists. It is given by 



(1 + A + A2).^1 - i(3 + A)r-i + ^(9 - 2A + X^)t-^ + ,^(3 + \^)t~^ - 1 - A + A^ 
= A[i(2 + 3A + A2 + A3)r-i-4A] " ^^^^ 

The critical density is 7]c = ri+{Tc), i.e. 

A[i(2 + 3A + A2 + A3)r-i-4A]" ^ ^ 

As A increases, Tc increases, while rjc decreases. Equations (^)-(^) generalize to A > 1 the results predicted by the 
PY solution for SHS,lI 



V9-6r-i+r-V2 + r-i-3 ^, „_3 



7r-i - 12 



= 3(2 + V2) ~ 10.24, ?7c = -V2- 2 ~ 0.1213. (35) 



The lack of a vapor branch of the spinodal line is a feature that our model inherits from the PY-SHS solution. ^3 
Our model also has in common with the SHS limit, as well as with the solution of the PY equation for finite A,Ej 
the existence of regions in the temperature-density plane, inside which the physical quantities cease to take real 
values. According to Eq. (p7|), this happens when 1 + Pirj + /32??^ + /^s??"^ + f^irj'^ < 0. Let us call T]i{T) (i = 1, . . . , 4) 
the four roots of the quartic equation 1 + fiirj + /92'7^ + /^s??^ + PiV"^ = 0, with the convention that, whenever they 
are real, ryi < 772 < '73 < ?74- It turns out that the roots 771 (T) and '72(1^) are real only if r is smaller than a 
certain threshold value rth (or, equivalently, if T < Tth); they define a dome-shaped curve with an apex at a density 
^th = ?7i(21:h) = ?/2(21:,h). Analogously, the other two roots, risiT) and rii{T), are real only if r < t^^ (T < T^^) and 
define another dome-shaped curve with an apex at r/^j^ = 'ysl^th) ~ '?4(^th)- Consequently, the parameter L2 becomes 
complex inside the intervals r?i(T) < 77 < r72(T) (for T < Tth) and r?3(T) < 77 < r74(T) (for T < T[^^). However, the 
existence of the second region is a mathematical artifact since it affects unphysically high densities. For instance, for 
A = 1.125 we have T^^ = 0.535147 and 7]'^^^ = 0.78692, this density being larger than the one corresponding to the close 
packing value rjcp = V^tt/G ~ 0.740. In fact, in the SHS limit (A 1) the second region collapses into the line 77 = 1 
and disappears. On the other hand, the region 771 (T) < rj < 772 (T) is inside the curve 77 — rj±{T) and persists in the 
SHS limit. In the latter case, the threshold values coincide with the critical ones, i.e. Tth = ''"c: %h = Ic- For finite A, 
Tth is smaller but practically indistinguishable from Tc and ryth is slightly smaller than 77c. In the illustrative case of 
A = 1.125, we have Tc = 0.502324, 77^ = 0.11299, Tth = 0.502269, and 7?th = 0.10980. The existence line, p = pi,2(r), 
and the liquid branch of the spinodal line, p = p^{T), are shown in Fig. |^ for A = 1.125. 
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FIG. 7. Liquid branch of the spinodal hne (sohd hne) and existence hne (dashed line) for A — 1.125, according to the present 
model. The model does not give real values for the physical quantities inside the grey region enclosed by the existence line. 

Recently, it has been suggestedll that, while the critical point is in general very sensitive to the range of the 
interaction, the second virial coefficient has a fairly .canstant value at the critical temperature, i?2(?c) ~ ^t^- In the 
case of the SW interaction, computer simulationgc3C3 show that the value of the second virial cpefficient B2{T) = 
— 27r[a::(A"^ ~ 1) ~ l]/3 at T = Tc is much less sensitive to the width A than the critical temperature.Ej^Ej Moreover, the 
PY-SHS solution predicts B2{Tc) = -(2 + 3\/2)7r/6 ~ -1.047r. Thus, by assuming that B2iTc) is independent of A 
and is equal to its value in the SHS limit, the criterion of Ref. ^ allows one to estimate the critical temperature as 

T, = l/ln[l + 3(2 + %/2)/4(A3-l)]. (36) 

Since the model presented in this paper is constructed as a simple generalization of Baxter's solution of the PY 
equation for SHS, we expect its predictions for the critical values of the temperature (Tc), the density {pc), and the 
compressibility factor (Zc) to be more reliable for small values of A — 1 than for wide wells. Ijtowever, to the best 
of our knowledge, the simulation estimates for those quantities are only available for A > 1.25,E3 and so we cannot 
make a comparison for smaller widths. Table || shows the values of T^, Pc, and Z,, for several valu es i pf iA, as estimated 
from computer siiiLulations,E3 and as predicted by the model, by Eq. (|36|), by perturbation theory,Ejl23 and by the PY 
integral equationO 



TABLE I. Critical constants of the square-well fluid for several values of the width parameter A. 


A 








Source 








1 


Qa 


0.232 


0.379 


PY equation^ 






1.1 


0.455 
0.461 


0.219 


0.372 


This wo 
Eq. ^) 


s 


1.125 


0.502 
0.512 
0.594 
0.587 


0.216 

0.46 
0.71 


0.370 

0.42 
0.74 


This wopk 
Eq. 

Perturbation theory 
Perturbation theory 


i3 




1.25 


0.764 
0.729 
0.766 


0.370 
0.203 


0.29 
0.360 


Comput 
This wo 
Eq. M) 


er simulation^ 



12 



0.913 
0.850 



0.34 
0.48 



0.43 
0.47 



Perturbation theory 
Perturbation theory 



1.375 



0.974 

0.960 

1.046 

1.11 

1.08 



0.355 
0.193 

0.34 
0.36 



0.30 
0.349 

0.39 
0.40 



Computer simulatio 
This woris 

Eq. 

Perturbation theory 
Perturbation theory 



s 



1.5 



1.219 

1.209 

1.367 

1.205 

1.35 

1.33 



0.299 
0.184 

0.200 

0.31 

0.29 



0.30 
0.339 

0.37 
0.36 
0.37 



Computer simulatio 
This wock 

PY equationO j-. 
Perturbation theory 
Perturbation theoryEi 



1.625 



1.479 
1.738 
1.70 
1.61 



0.177 

0.27 
0.26 



0.330 

0.38 
0.36 



This wopk 
Eq. @y 

Perturbation theory 
Perturbation theory 



a 



1.75 



1.811 

1.777 

2.164 

2.04 

1.93 



0.284 
0.170 

0.25 
0.24 



0.35 
0.322 

0.38 
0.36 



Computer simulatio 
This wopk 
Eq. 

Perturbation theory 
Perturbation theoryEi 



1.85 



2.036 
2.550 
2.33 
2.23 



0.165 

0.25 
0.23 



0.317 

0.37 
0.35 



This wopk 
Eq. ||)y 

Perturbation theory 
Perturbation theory 



2.764 

2.466 

3.208 

2.88 

2.79 



0.225 
0.159 

0.24 
0.23 



0.32 
0.310 

0.37 
0.35 



Computer simulatio 
This wof44 
Eq. ©y 
Perturbation theory 
Perturbation theoryEi 



= 0.0976 
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Comparison with MC simulation data shows that the critical temperature predicted-by our model, Eq. ( |32D , tends 
to be smaller thsmptiie correct value, while Vliegenthart and Lekkerkerker's criterionJ£3 as well as second-order per- 
turbation theoryHa tend to overestimate it. What is indeed remarkable is the fact that Eq. ( ^2| ) provides the best 
agreement with computer simulations for A < 1.75 (except at A = 1.25, in which case Eq. (^) is better). In the case of 
the critical value of the compressibility factor, the best agreement corresponds to our model, except at A = 1.75, where 
Chang and Sandler's perturbation theory gives a better result. On the other hand, the critical densities predicted 
by our model are around 30-45% smaller than the simulation values, a discrepancy that can be traced back to the 
solution of the PY equation for SHS. Since the computer simulations and all the theories share the property that pc is 
a monotonically decreasing function of A, we can conclude that the correct value of pc in the limit A — s- 1 is certainly 
larger than the simulation value pc = 0.370 at A = 1.25, while Baxter's solution yields pc = 0.232. Thus, we can 
expect this failure to reproduce accurately the critical density is also present in the PY|-approximation for finite A. 
This is confirmed by the results of a numerical solution of the PY equation for A = 1.5,E3 which gives a value of pc 
rather close to the one obtained here. 

Now let us compare the general density dependence of the compressibility factor predicted by the model with 
available computer simulationSr— .We start with the smallest value of A that, to our knowledge, has been analyzed in 
simulations, namely A = 1.125.E3 Figure ^ shows Z{p) for A = 1.125 and T = 0.5, 0.67, and 1. Strictly speaking, the 
curve representing the model at the lowest temperature corresponds to T = Tc ~ 0.502 rather than to T = 0.5. This 
is because at T = 0.5 there exists a small density interval around pc — 0.22 [cf. Fig. ^ where the model gives complpj. 
values. The next width we consider is A = 1.4, for which extensive molecular dynamics simulations are available.Ej 
The results for T = 1.25, 1.43, 2, and 5 are plotted in Fig. ^. We observe that, except at the highest temperatures 
(T = 1 for A = 1.125, T > 2 for A = 1.4), our model presents a better general agreement with the simulation data 
than the TL perturbation theory. 
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FIG. 8. Density depe»4ence of the compressibility factor for A = 1.125 and three temperatures. The circles are the results 
of computeii-simulations,Ej the solid lines are the results from our model, and the dashed lines are the TL perturbation theory 
predictions. lZI 
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FIG. 9. Density dependence of the compressibility factor for A = 1.4 and four temperatures. The circles correspond to an 
empirical formula fitted to molecular dynamics rssMlts,L3 the solid lines are the results from the present model, and the dashed 
lines are the TL perturbation theory predictionsEZI 



IV. CONCLUDING REMARKS 



The main objective of this paper has been to propose an analytical model that could be useful to describe the 
structural and thermodynamic properties of systems, such as colloidal dispersions, composed of particles effectively 
interacting through a hard-core potential with a short-range attraction. As the simplest interaction that captures 
both features we have considered the square- well (SW) potential (|l|). On the one hand, perturbation theory becomes 
unreliable when the range A — 1 of the attraction is small, especially at low temperatures, i.e. when a certain degree 
of "stickyness" among the particles becomes important. On the other hand, the widely used sticky-hard-sphere 
(SHS) interaction model combines temperature and well width in one parameter only, thus lacking the flexibility to 
accommodate additional changes in width and/or temperature. Our approach intends to fill the gap between these 
two theories. 

The model presented in this paper is based on the one proposed in Rcf. where the functions R(t) and R{t) defined 
by Eqs. (j^) and ( ]l5| ) were approximated by rational forms, Eq. The parameters in these functions are constrained 
to yield a finite value for the isothermal compressibility by Eqs. (|19D- ([2^). This still leaves two free parameters, A and 
T, the latter being defined in Eq. (p5|), as unknown functions of the packing fraction n, the temperature parameter 
X = e^^^ — 1, and the well width A. Apart from their zero-density limits ( |2^ ) and (p6|), two extra conditions are 
needed to fix those parameters and close the construction of the model. Thus, depending on the physical situation 
one is interested in and/or on the degree of simplicity one wants to keep in the model, it is possible to choose different 
criteria to determine A and r. For instance, one could impose cectain continuity conditions on the cavity function 
y(r) = g{r)e^^'^^/^^"'^ at the points where the potential is singularj^ alternatively, one could require thermodynamic 
consistency among the virial, compressibility, and energy routes. Of course, other choices are possible. In the original 
formulation of the model, A was assumed to be independent of density (hence A = x) and the second condition was 
the continuity of y(r) at r = Act, this giving rise to a transcendent equation that needed to be solved numerically. On 
the other hand, in this paper we have simply assumed that both A and r are independent of density, Eq. (28). The 
assumption for r is expected to be especially adequate for narrow potentials, since in the SHS limit the role of r is 
played by the parameter tshs = [12a;(A — 1)]~^, which is indeed independent of density. With these choices for A and 
T the problem remains fully algebraic and all the parameters can be expressed in terms q£ the solution of a quadratic 
equation, in analogy with what happens in Baxter's solution of the PY equation for SHSJil In fact, the model includes 
such a solution as a limit case. Given the scarcity of simulation results for narrow wells, we have been forced to carry 
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out a comparison for cases with A > 1.125. In spite of that, the results show a very good general performance of 
the structural and thermodynamic properties predicted by the model, correcting the inadequacy of the perturbation 
theory predictions in the low-temperature domain. It is interesting to remark that the model provides an explicit 
expression of the critical temperature as a function of the well width, Eq. (|3^), which is accurate even for rather wide 
wells. 

The results of this paper could also be useful in connection with a recently proposed extension of the law of corrje^ 
spending states for systems, such as colloidal suspensions, that have widely different range of attractive interactions.ESl 
Given an interaction potential = (^i.op(r) + ipiittif), where (prcp and (^att are the repulsive (not necessarily hard- 
core) and attractive parts, respectively, one can define a (temperature-dependent) effective hard-core diameter 



an effective well depth 



cr = / dr 



e = -¥?(ro). 



1-e 



d(f{r) 



dr 



(37) 



(38) 



and a (temperature-dependent) effective well width 



A = 



1 + {B; - 1) (l - e'^^^^ 



1/3 



(39) 



where 



dr r 



-v{r)/kBT 



(40) 



is the reduced second virial coefficient. Then, according to the extended law of corresponding states,E3 the com- 
pressibility factor for a wide range of colloidal materials is a function of only the reduced temperature T* — kBT/t, 
the reduced density p* = pa^, and the range parameter A, i.e. Z = F{T* , p* , X), where the function F is hardly 
sensitive to the details of the potential. As a consequence, an accurate prescription for the function F based on 
the SW interaction for variable width can be used to determine the thermodynamic properties of a wide class of 
colloidal suspensions. As an illustrative example. Fig. ^ shows the dependence of the reduced critical temperature 
T* = kBTc/e on the effective width A for several interaction potcntials.E3 The prediction of our model, Eqs. ( p9| ) and 
(P^), is also plotted. The extended law of corresponding states work very well up to A ^ 1.3. For larger interaction 
ranges the values of T* for the generalized Lennard- Jones and the hard-core Yukawa potentials tend to lie slightly 
above those corresponding to the SW interaction. 



16 



2.0 




0.0 u ' ^ ' ^ ' ^ ' — I 

1.0 1.2 1.4 1.6 

X 

FIG. 10. Dependence of the reduced critical temperature on the efective well width, as obtained from|-Gpiwjiter sim- 
ulations for several interaction potentials: square- well (solid circles) ,l3 hasjd-core Yukawa (open squares) jHfl'Ej'H general- 
ized Lennaj[(i| Jones (open circles) ,Ej'E3 a-Lennard-Jones (open diamond) ,E3'lj an effective colloid-colloid intep|q ti| i pn. (up 
triangles) J2£ICj and an effective interaction for non-additive mixtures of asymmetric hard spheres (down triangles) .Ej'l3H The 
solid line is the prediction of our model, Eqs. (^9]) and (^^. 

In the near future, we plan to extend the model presented in this paper in two directions. First, the case of a 
mixture of particles interacting through SW potentials with different values of ct, e, aud A will be analyzed. This 
generalization must be such that one recovers the cases of a mixture of hard spheres ,EJ a mixture of sticky hard 
spheres,E3 and a monodisperse SW system in the appropriate limits. As a second extension, we will study the case 
of an interaction model made of a hard core plus a square shoulder plus a square well, for which, in_addition to the 
conventional gas-liquid phase transition, a liquid-liquid transition in the supercooled phase appears.EEl Ej 
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APPENDIX A: SOME EXPLICIT EXPRESSIONS 
1. Expression for ^(O) 

In this Appendix we list some of the expressions that are derived from the model. The long-wavelength value of 
the structure factor is obtained by using the model ( p^ ) in the general expression ( p^ ) . The result is 

S{0) = 5(YT2^ " ^""^ [1 + (2 + A)^' - 3(1 + X)L2] 

+2ri^ [15 - 6(14 - A + 19X^ - A^ - X^)L2 + 120(1 + A + X^)LI 
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+ (50 + 16A + 27A2 - 2A3 - A'*) A' - 30(5 + 4A + ?,\^)A'L2 
+15(3 + 2A + A2)A'^ 

-277^ [10 + 3(7 - 53A - 13A2 + 7A^ - 8A'')i2 

-60(1 + A)(l - 4A + A2)(1 + a + \^)Ll 

+ (19 + 59A + 9A2 - A^ + 4A^) A' 

+3(11 - 67A - 114A2 - 66A3 - 13A4 + 9A^)^'L2 

+3(3 + 2A + + 7A + 3A2 - A^)^'^ 

+'n'^ [5 - 12(1 + A + A^ + llA^ - 4A'')L2 + 240A'''(1 + A + \^)LI 

+2 (7 + 8A + 9A2 + lOA'"^ - 4A'') A' - 12(1 + 3A + GA^ + 24A''' + 17A'* + 9A^)A'L2 

+3(3 + 2A + A2)(1 + 2A + 3A2 + AX^)A''^ | , 



where we have made use of Eqs. (19)-(E' 



2. Expressions for the coefficients in Eq. (27) 



Equation (E5[) reduces to a quadratic equation for L2. Its physical sohition is given by Eq. (27), where 



ai = 2A'(2 + A) + i (1 + 4A + A^ - 3A') r-\ 

D 

a2 = 3 + ^'(7 + 2A - 3A2) - Y2 + A + IBA^ + A'(23 + 15A + 15A2 + 7A^)] 

as = -2 - 2A'(1 + 2A + iX^) + -\7 + \-2\^ + A'{7 + 15A + 21A2 + llA^ + 6A'')] t'^ 

6 

/3i ^ -4 - 4A'(2 + A) + i(5 + 2A - A^ + ?,A')t-^ - \t-^ , 



/32 = 6 + 6A'(3 + 2A + A^) + 4A'^(2 + A)^ 

^ [9(3 + A) + A'(59 + 51A + 3A2 - 5A^) + 12A'^(2 + A) 



6 

1 r 
36 



31 + 8A + 18A2 - 4A^ + A^ + 12A'(5 + A) + QA'" 



/33 = -4 - 12A'(1 + A + A^) - 4A'^(2 + A)(l + 2A + 

+ 1 [3(4 + A + A^) + A'(29 + 36A + 30A2 + - ?,X^) 

+A'^{17 + 43A + 30A2 + A''' - A'') 

[29 + 13A + 27A2 + 7A-'^ - AX'^ + yl'(68 + 80A + 80A2 + 16A-'^ + 7A^ + A^) 
36 '■ 

+3A'^(13 + 13A + A2-3A3) 



/34= [l + A'(l + 2A + 3A2)]'-i [7 + A + 4A2 

+2A'(7 + 15A + 15A2 + 5A=^ + eA*) + A'^{1 + 2A + 3A2)(7 + 15A + SX"^ - A^) 
1 



144 



[49 + 38A + 9A2 + 8A-'^ + 16A* + 2A'(49 + 88A + 112A^ + 80A^ + 35A'' - 4A^) 



-A'" (49 + 138A + 219A^ + 124A^ + 27A'' + 18A'^ + A^) r^^ 



18 



70 = 1 + A- \ 
6 



71 =2(l + A-A2)--(3 + A-^)r-\ 



1 



-4A + ^(2 + 3A + A^ + X^)t~^ 



72 = A 

The other root is incompatible with (E3[) and then must be discarded 



(A9) 
(AlO) 
(All) 



APPENDIX B: LOW-DENSITY BEHAVIOR OF THE MODEL 



In the low-density limit, Eqs. (|l9|)-(22) become, respectively, 



Si = 0{rj), 



(Bl) 
(B2) 



S2 = ~-[l-xiX'-l)]+0{7^), 



53 = -^ + i[l-x(A^-l)]+0(77). 



To first order in density, Eq. (]27|) gives 



where 



l'^P = --xX{\ - 1) {l + a;(A - 1)^ [i _ 2A - A^ - 4a;A(l + A)] } 



Substitution into Eq. ( p^ ) yields, after some algebra, 

Co + Ci(A-l) + C2(A-l)2" 



2C2 Ci+2C2(A-1) 



i3 



where we have called 



C2 = -3a;(l + a;)(A2 - 1), 



Ci = 2x{l + x){X- 1)2(1 + 2A), 



(B3) 
(B4) 

(B5) 
(B6) 



(B7) 

(B8) 
(B9) 



Co 



(1) 



A- 1 2 



x{2 - SA^ + 3A'') - 2a;'(A - 1)^(2 + 4A + 3A^) - 3x-'(A^ - 1) 



(BIO) 



The expression for the exact radial distribution function, gcxact('')i to first order in density was derived by Barker and 
Henderson for the case A < 2.CJ The corresponding expression for inexact (i) can be found in Ref. ^ From Eq. (B7) 
one easily gets 
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g{r) - 5cxact(r) = -[C2{r - if + Ci(r - 1) + Co][e(r - 1) - Q{r - A)]r; + 0{if). 



(Bll) 



Thus, the difference (to first order in density) is non-zero in the interval 1 > r > A only. In particular, 

5(l + )-ffcxact(l+)=Co77 + 0(??'), 

where0 



(B12) 



5oxact(l ) = 1 + .T + 



^ + |(15 _ i6A3 + 6A'*) - 2x^{\ - 1)(4 + 4A + A2 - SA^) 



+3x3(A'- 1)2] 7^ + 0(772). 



(B13) 



Note that the relative coefficient Co/(l + a;) vanishes in the SHS limit . As said in Sec. ||, the model presented in this 
paper does not enforce the verification of Eq. (p4). In fact, Eq. (Bll) implies that 



5(A-) - (1 + x)g{\+) = A-i [C2(A - 1)^ + Ci(A - 1) + Co] 77 + 0{rf) 

= ^^4v^ [A(ll - A - A2) - x{X - 1)(A + 2)(3 + lOA - SA^) 

+6x2(A - 1)2(1 + 3A + 2A2)] rj + 0{rf). 

Finally, from Eq. (Al) or, equivalently, inserting Eq. (Bll) into Eq. (^, we get 

5(0) - S'cxact(O) = 2x(A - if [7 + 23A + 30A2 - - - x{\ ~ I) 

X (25 + 84A + 78A2 + - QX^) + l^x^{\^ - lf{l + 2A)] 77^ + 0(^3^^ 

where 

5'oxact(0) 1 - 8 [1 - x{\^ - 1)] 77 + 2 [17 - a;(A - 1)(19 + 19A + lOA^ + SIA^ - 3A'' - ZX^) 
-2a;2(A - 1)^(8 + 16A - SA^ - 38A3 - lOA'^) + l%x^{\^ - if] 77^ + 0(773). 



(B14) 



(B15) 



(B16) 



APPENDIX C: THE HARD-SPHERE AND STICKY-HARD-SPHERE LIMITS 



1. Hard spheres 

The SW potential becomes equivalent to the HS potential if A = 1 at any no n-zer o temperature T or if T ^ cxd at 
any width A. Let us first consider the latter limit. Making .t = in Eqs. ( A2)-( All ), one gets 



ai =0, a2 = 3, as = —2, 



(CI) 



As a consequence, Eq. 



with 



/3i=-4, /32 = 6, /33 = -4, /34 = 1, 

70 = 1 + A, 71 = 2(1 + A- A2), 73 = -4A2. 
|) implies that ^2 = 0. Thus Eq. ( p^ becomes 
1 1 + Lit 



Urj 1 + Sit + S2t^ + 53*3 ' 



(C2) 
(C3) 

(C4) 



Li 



1 + 27;' 



o 3 77 



21 + 277' 



(C5) 
(C6) 
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So = -- 



1 1-77 



2 l + 2r; 
(1 - ri f 



(C7) 



(C8) 



1277(1 + 277) ■ 

This is precisely the form adopted by F{t) in the analytical solution of the PY equation for hard spheresESElS The 
expression for 5'(0), Eq. (Al), simply reduces to 



S{0) 



(1 + 277)2- 

The case A = 1 is not considered in this Subsection, as it is a particular case of the SHS limit. 



(C9) 





= r-i 


a2=3-2T-\ 


Pi = 


-4 + 


2T-i-ir-2, (32 


/33 = 


-4 + 


6t-1-2t-2, (3i 




1 




70 = 2 


T 

6 





as = -2 + 2r-\ 



/34 = 1 - 2t 



73 



-1 i_ S -2 
6 ' ' 



-4 + -t" 



Therefore, Eq. ( p7| ) reduces to 

_ 1 - 77 (1 + 27,) [(1 - 77)2 + 27,(1 - 77)t-i - 177(2 - 577)r-2] 



-(l-77)(l + 277 - 7,r-i) 



2477 (l-77)(l + 277)- j3j(l + 47, - 14772)t-i 

Taking the limit A ^ 1 in Eq. (|l^) we get 



2. Sticky hard spheres 

Let us now take the limit x — > 00, A — s- 1, with x{X — 1) = fini te in t he m odel proposed in this paper. In that limit 
the parameter r is finite, cf. Eq. (28), while A' = 0. Equations (A2)-( A11 ) become 

(CIO) 
(Cll) 

(C12) 
(C13) 

(C14) 

(C15) 
(C16) 
(C17) 
(CIS) 

From Eq. (^) we have 
(C19) 



with 



m = - 

Li = 



1 l + Lit + L2t^ 
12?7 1 + Sit + S2t^ + S3P 

1 + 577 677 



1 + 2?7 1 + 277 



L2, 



3 V 



677 



2 1 + 277 1 + 277 

1 1 - 77 1 - 477 . 

2 1 + 27/ ^ 1 + 277'' 



(1-77)2 



1-77 



L2. 



1277(1 + 277) 1 + 27/ 

This coincides with the analytical solution of the PY equation for sticky hard spheres 

(1-77)2 



3(0) = 



[l- 77 + 1277L2]^ 



(1 + 277)2 

Of course, the results for hard spheres are recovered in the high-temperature limit (t^^ 



0). 
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